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SYMMETRIES OF 2-LATTICES AND SECOND ORDER 
ACCURACY OF THE CAUCHY-BORN MODEL 

B. VAN KOTEN AND C. ORTNER 

Abstract. We show that the Cauchy-Born model of a single- 
species 2-lattice is second order if the atomistic and continuum 

£N| kinematics are connected in a novel way. Our proof uses a gener- 

alization to 2-lattices of the point symmetry of Bravais lattices. 

Moreover, by identifying similar symmetries in multi-species 
pair interaction models, we construct a new stored energy density, 
3 using shift-gradients but not strain gradients, that is also second 

^H order accurate. 

^H These results can be used to develop highly accurate continuum 

£" — models and atomistic/continuum coupling methods for materials 

such as graphene, hep metals, and shape memory alloys. 

< 

c3 The Cauchy-Born model is a widely used continuum model for crys- 

tal elasticity |4[|5J|20] . Moreover, it is a crucial ingredient in a new class 
of atomistic/continuum multi-scale methods (3, 12} 13}[T9j[20). Formal 



^1" 



a 



1. Introduction 



considerations and rigorous analyses have shown that the Cauchy- 
Born model for Bravais lattices (simple lattices) is second order ac- 
curate [4j[8J[lO}[l6] . By contrast, one expects that its generalization to 
multi-lattices should be only first order accurate (8) , due to the absence 
of point symmetry in general multi-lattices. 

In the present work, we identify two non-trivial generalizations of the 
Bravais lattice point symmetry. We show that this leads to second order 
accuracy of the classical Cauchy-Born model for a single-species 2- 
>• lattice, provided the atomistic and continuum kinematics are connected 

in a novel way. Moreover, we identify a new stored energy density for 
general multi-lattices under pair interaction model, and we show that 
this energy is also second order accurate. 

While these are interesting observations on their own, they have im- 
portant consequences for computational materials modeling. For exam- 
ple, our results provide higher-order continuum approximations without 
requiring the use of ^-conforming numerical methods. (Higher order 
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methods" . 

1 



2 B. VAN KOTEN AND C. ORTNER 

continuum models can always be constructed provided they are dis- 
cretized using higher-order conforming numerical methods (Tl.) Our 
own interest in this issue is the application of the Cauchy-Born model 
in atomistic/continuum coupling methods. At the atomistic/continuum 
interface, the order of accuracy is typically lower than in the continuum 
bulk, due to the loss of interaction symmetry [7,15,18,22]. By employ- 



ing a second order continuum approximation, the loss can be made less 



severe. For example, in the blended quasicontinuum method 22 , the 
interfacial error is controlled in terms of an approximation parameter 
k called the blending width. In ID, it is shown in [22] that the leading 
order term in the error decreases as k~ 2 when the site energy has point 



symmetry (see Definition (15) below). By contrast, we show in 17 
that if the site energy does not have point symmetry, then the error 
decreases as k~ a. In higher dimensions, the loss in accuracy can be so 
severe that the method may cease to be consistent [17]. In addition, 
the symmetries that we identify allow us to investigate ghost-force re- 



moval techniques as discussed in 18, 19 , which cannot be employed 



when site energies do not possess point symmetry. 

Our framework encompasses all single-species materials with 2-lattice 
structure and general multi-lattices modeled by pure pair interactions. 
Physical systems of interest that can be described within this frame- 
work are hep metals (e.g., Mg, Ti, Zn), the honeycomb lattice (graphene) 
or the diamond cubic structures (e.g., diamond or silicon). Common 
multi-species materials with multi-lattice structure are shape memory 
alloys such as Ni-Ti, Fe-Ni, Ni-Al, which are often modelled using only 
pair interactions [9,11 . 



Outline. In Section[2j we introduce the atomistic energies that we con- 
sider, and the associated notation for atomistic kinematics. In Section 
|3j we derive the classical Cauchy-Born energy for simple and multi- 
lattices and prove first order accuracy of the energy in the multi-lattice 
case and second order accuracy in the simple lattice case. The second 
order result is based on the point symmetry of Bravais lattices, which 
motivates the search for similar symmetries in multi-lattices in Section 
[1] We identify two non-trivial types of point symmetries in this sec- 
tion. In Section [5j we exploit the symmetry in single-species 2-lattices 
to prove second order accuracy of the classical Cauchy-Born energy. 
This result does not extend immediately to the case of multi-lattice 
pair interactions. Instead, in Section [6] we construct a novel stored 
energy density for which we are again able to exploit the symmetries 
discovered in Section [4] to prove second order accuracy. 

For the sake of simplicity, we formulate all our results for 2-lattices. 
However, it is straightforward to see that the results for pair interac- 
tions (and only these) extend to general multi-lattices. 
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2. Atomistic models for 2-lattices 

2.1. Atomistic kinematics. For d G {1,2,3}, a rf-dimensional Bra- 
vais lattice (simple lattice or 1-lattice) is a set of the form BZ d for some 
strain B G GL(cf). A d- dimensional 2-lattice is a set of the form 

(BZ d + Po )U(BZ d + Pl ), (1) 

for some shifts p = {po,pi} G M d x M. d . 

We call Z d the reference lattice and a point (GZ^a szte. A conve- 
nient index set for a 2-lattice is 

A :=Z d x {0,1}, 

which we call the reference list. The reference list serves a purpose 
similar to the reference domain in elasticity. We think of the elements 
of the reference list as atoms or nuclei, and we call (£, a) the atom of 
index a at site £. 

An atomistic deformation is a map y : A — y M> d . We will use the 
notation y a (£) :— y(£, a) for evaluations of functions defined on A, and 
we call y a (0 the deformed position of atom (£, a). To make our analysis 
as simple as possible, we impose periodic boundary conditions on the set 
of deformations. Fix iVeN. We call a map u : Z d — y M. d an N -periodic 
displacement if 

«*(£) =«a(^ + ^) for all f, r? G Z d , aG {0,1}, 

and we call y an N -periodic deformation if for some A^-periodic dis- 
placement u and some strain B G GL(d) we have 

y a (0 = B£ + u a (0. 

Throughout the remainder of the paper, we will assume that all defor- 
mations y are A^-periodic. 

Remark 2.1 (Index versus species). It is important for our purposes 
to distinguish between lattices that are composed of identical atoms 
and lattices that are composed of atoms of two or more species. Thus, 
we will draw a careful distinction between the species of an atom and 
its index. The species is the type of atom, e.g. Cu, Zn, C. The index 
belongs to the set {0, 1}, and it tells us which of the component Bravais 
lattices making up the 2-lattice should be associated with the atom. 
We assume that atoms of the same index must be of the same species. 

2.2. Atomistic energies. Let y be an A^-periodic deformation, and 
let f2jv := {0, 1, . . . ,N — l} d be a periodic cell. For a pair interaction 
model, the atomistic energy takes the form 

^(i/) : =E|E E IMIwfa)-i/«(fll) 

££Cl N ae{0,l} (r;,/3)eA 
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where the functions 4> a/ 3 : (0, oo) — > M. are pair potentials that may de- 
pend on the species of interacting atoms. If the material is composed 
of atoms of a single species, then there should be no dependence of the 
potentials on a and 0. Examples of multi-lattice pair interaction mod- 
els describing interesting mechanics such as the shape-memory effect 
are described in (9,11 



We call the inner sum 

E E sMlvM - v*®\) (2) 

ae{o,i} (r),p)eA 

the site-energy at £. To simplify our analysis, we will assume that 
(pa/3 G C 3 ([0, oo); E), even though physically realistic potentials (such 
as the Lennard- Jones potential) may have singularities at 0. We also 
assume that for some cut-off radius r c e (0, oo), 

0«/3( r ) = whenever r > r c . 

Thus, the sum defining the site-energy is finite as long as y(A) does not 
have an accumulation point. To avoid discussing this purely technical 
point, we assume in the following that the interaction range is finite in 
the reference configuration. 

In our analysis, we will also allow a more general class of potentials 
than pair interactions. We only require that the total energy can be 
decomposed into a sum of localised site-energies. To that end, let 

K C Z d x {0, 1} x {0, 1} \ {(0; 0, 0), (0; 1, 1)} 

be a finite interaction range. Throughout the remainder of the pa- 
per, p = (p; a, (3) will denote an element of 1Z. Given p e 1Z and a 
deformation y e &, we define the 2-lattice finite difference 

and the ft-tuple D n y{£) := (D p y(£,)) pen . 

We assume that the total atomistic energy takes the form 

^(y) = E v (Dny(0), (3) 

£en N 

where V : (l^) 7 ^ — > R is a site potential. In our analysis, we will assume 
that V E C 3 ((M d ) n ;'K). Clearly, the pair interaction site-energy ([2| is 
of this form. 

Moreover, typical EAM potentials for hep metals |2| or bond-angle 
and bond-order potentials for carbon structures (graphene) (6,21 



can 



be written in this form. However, it may be impossible to express po- 
tentials arising directly from quantum mechanics or electronic structure 
models as in pi). 

Remark 2.2. Fixing the interaction range 1Z in the reference domain, is 
justified by the fact that we consider only elastic effects in this paper. 
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3. The Cauchy-Born energy 

The Cauchy-Born energy is an elastic energy which provides a good 
approximation of (|3~| for deformations which are close to homogeneous 



(i.e., smooth) EM, 16 



3.1. Continuum kinematics. In continuum models for 2-lattices, the 
kinematic variables are a deformation field Y G C 1 (M d ; IR d ) and a shift 
field P G C°(R d ; R d ). We say that a pair of fields (Y, P) is iV-periodic 
if, for some macroscopic strain B G GL(cf), 

P(x + Nrj) = P(x) and y(x + Nrj) = NBr] + y(x) for all r) G Z d . 

3.2. The Cauchy-Born energy for Bravais lattices. We first re- 
view the Cauchy-Born approximation for Bravais lattices. In this case, 
we may ignore the shifts, and hence atomistic deformations are now 
maps from Z d to M d , and the interaction range 1Z is a subset of Z d . 
The continuum kinematic variable is just a single iV-periodic deforma- 
tion field Y G C 1 (R d ] R d ). 

Set Q := [0, l) d . We observe that Nfl is a periodic cell for an N- 
periodic continuum deformation. The Cauchy-Born energy (for Bra- 
vais lattices) takes the form 

g c (Y) = I W{VY) dx, (4) 

JNn 

where W : M> dxd — > R U {+00} is the Cauchy-Born strain energy den- 
sity. For F G GL(<f), W(F) is defined to be the atomistic energy per 
unit volume in the lattice F1, d . That is, for the atomistic deformation 

y F (0 ■■= F£, 

W(F) := lim N- d V V(D n y F (0) = V(D n y F (0)), (5) 

£en N 

since D n y F (0 = D n y F (0) = {Fp}^ for all f G Z d . 
Let Y be a continuum deformation. We observe that 



W(VY(x)) = V(V n Y(x)), and S C {Y) = / ^(V^Y) da;, (6) 

where V n Y(x) := {V p F(x)} p ^. 

3.3. The Cauchy-Born energy for 2-lattices. We now explain how 
the Cauchy-Born model is traditionally generalized to 2-lattices. Let 
Y and P be iV-periodic deformation and shift fields. The Cauchy-Born 
energy (for 2-lattices) takes the form 



S C (Y,P) : = / W(VY,P)dx, 

Jnq 

where W : M> dxd x (IR d )™ — > IRU{+oo} is the Cauchy-Born strain energy 
density for 2-lattices. 
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As in the case of Bravais lattices, W(F, P) is the energy per unit 
volume in a lattice subjected to a homogeneous deformation with strain 
F and shift P. That is, for the deformation y F)P defined by 2/„' p (£) : = 
Ft; + aP, we have 

W(F,P) := lira IT* V V(D n y F > p ($) = V (D n y F ' p (0)) , (7) 

since D n y F ' p (£) = D n y F > p (0) = {Fp + (J3 - a)P} pe7l for all £ G Z d . 
(Recall the convention p = (p; a, /?).) 

We now generalize (J6J). For p <E 1Z, we define the 2-lattice directional 
derivative V p by 

V P (Y, P)(x) := V P F (x) + (/3 - a)P(x), (8) 

and we set V n (Y,P)(x) := {V P (F, P)(x)} p£7r Thus, by @, 
W(VY, P) =^(V*(Y, P)(z)), and 
g\Y,P)=l V(V n (Y,P)(x))dx. (9) 

3.4. Error estimates for the Cauchy— Born energy. In order to 
compare the atomistic and Cauchy-Born energies, we must specify how 
to generate an atomistic deformation y from a continuum deformation 
field Y and shift field P. In this section, we adopt the classical ap- 
proach, however, we will see in Section [5] that for a different identi- 
fication of the atomistic and continuum variables, the Cauchy-Born 
energy can be a better approximation. 

The Cauchy-Born energy accurately approximates the atomistic en- 
ergy when the deformation and shift fields are "smooth" or "nearly 
homogenous." To make this precise, we define a family of increasingly 
smooth deformations Y N and shifts P . Let Y and P be fixed, 1- 
periodic deformation and shift fields. Define the scaled deformation 
and shift fields Y N and P N by 

-Nf\ ah/ ( X \ „„j tdN' _/# 



Y»{z):=NY{jj) and P» (x) := P (-) , (10) 

and define a corresponding atomistic deformation y N by 

y£(0--=Y N (0 and y?(0--=Y N (Z) + P N (0- (n) 

Observe that Y N , P N , and y N are all iV-periodic. We adopt the 
convention that the atomistic energy of an iV-periodic deformation is 
always the energy of the periodic cell Qjy, and that the Cauchy-Born 
energy is the energy of NQ\ that is, 

^(V N ) ■= E V ( D nV N (0), and 

(en N 

go(Y N ,P N ):= [ W(VY N ,P N )dx. 
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Since we have VF w (i) = VY (|r) , the strain energy density satisfies 
W(VY N (x),P N (x)) = W (vY (il) ,P (£)) , 

and so the Cauchy-Born energy has the scaling invariance 

N~ d g c (Y N ,P N ) = <f c (Y,P). 

This suggests that we should treat <f c (Y, P) as an approximation of the 
atomistic energy per atom N~ d S' aj [y N ^ . Indeed we show in Proposition 
that 

lim N- d ^{y N ) =£ C (Y,P). 

That is, <f c (Y, P) may be understood as the continuum limit of N~ d S' & (i/ N ^ 
as iV — > oo; see El. The proof is elementary, but we provide it never- 
theless for comparison with our subsequent analysis. 

Proposition 3.1 (Convergence of Cauchy-Born energy for 2-lattices). 
Let (Y,P) be a 1-periodic deformation and shift field. Let y N be the 
scaled atomistic deformation defined by (11). Then we have 

\ N - d ^(y N ) - <? C (Y,P)\ < CN- 1 {HWIU + IIVPIU} , 

where the constant C is a function oflZ and ||V||ci. 

Proof. Since the midpoint rule is exact for constant functions, a stan- 
dard quadrature estimate gives 

N -dgc(Y N ,P N ) = N- d [ V(V n {Y N ,P N )) dx 
= N- d J2v(Vn(Y\P N ))(0 

+ 0(\\V\\ c i{\\V 2 Y N \\ +\\VP N \\ }). 

\H ■ ■ * — l- 1 1 lloo II llooJ/ 

Thus, 

Err(N) :=N- d {g a (y N ) -g c (Y N ,P N )} 

=N~ d £ V(D n y N (t))-V(V n (Y N ,P N ))(0, 

+ 0(\\V\\ c i{\\V 2 Y N \\ +\\VP N \\ }). (12) 

\>> ' ' *— ' L I I lloO II MOO J / V ' 

By the mean value theorem, 

\V(D n y N (0)-V(V n (Y N ,P N )(0)\ 

<Y,\\V\\ci\D p y N (Z) -V P (Y N ,P N )(®\, (13) 
pell 
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and using Taylor's theorem, 

D P y N (0 - V P (Y\ P N ) (0 = ${£ + p)- y^(0 

- V p Y N (0 + (a- f3)P N ($ 
= Y(Z + p)-Y(0-V p Y N (0 

+ P{P N (£ + p)- !>"(£)} 
= 0(||V 2 F iV || + ||VP Af || ). (14) 

Ml II OO : ! II OO/ V ' 



Combining (12), (13), and (14) shows 



N- d \SHy N )-S c (Y N 1 P N )\<C{\\V 2 Y N \\ + \\VP N \\ } 

1 v*7 / \ /1 — L|| Uqq 11 I loo J 

where the constant C is a function of ||V||cri and 1Z. We now observe 

72\rN\\ _ at— l|IV72vll „„J 1 1 V7 D-WI I _ AT— 1 1 



V z y JV =rWTL and VP JV = iV" 1 VP 



OG 



OO) 



and the result follows. □ 



For Bravais lattices, the estimate given in Proposition |3. 1 1 can be im- 
proved if we assume that the site potential V has point symmetry [To] : 

V(g) = V(-{g- p } pen ) for all g G (K d )*. (15) 

We show in Lemma |3.1| that under physically reasonable assumptions, 
one can always take V to be point symmetric. Let R : M. d — > M. d defined 
by Rx = —x be point reflection in the origin. The first assumption is 
that 

£*(y) = £*( y o R) for all ?/6f. 

Physically, this assumption is motivated by the observation that per- 
mutation or relabelling of the atoms does not change the energy, and 
that Bravais lattices are invariant under the point inversion R, which 
means that R provides such a relabelling. 
The second assumption is that 

S & {y) = g*(-y) for all ?/ef, 

which is motivated by the principle that the energy should be un- 
changed if the configuration of the atoms is reflected. 

Remark 3.1. Permutation invariance can be expected to hold whenever 
there is only one species of atoms. If there is more than one species, 
then we can only assume that the energy is invariant under those per- 



mutations which preserve the species of the atoms. Hence, Lemma [3TT 
cannot be immediately applied to general multi-lattices. 

Lemma 3.1. Assume that for all deformations y, 

£*(y) = £*(-yoR), (16) 

where Rx = —x; then 
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where TZ = TZU -K, and V : (M^) 71 ->> R rfe/wied fry 

^({<?pW) := \V{{g P } P en) + \V({-g- p } p ^), 
is point symmetric (115]). 



Proof. First, we note that, if y is an iV-periodic deformation then so is 
— y o _R, since for all r\ G Z d 

-3/ o i?(£ + Afy) = -y(-£ - N V ) = -y(-0 = -y o R(£). 



By (16), we have 

^{y) = \^{y) + \^{-yoR) 

= \ E v ( D -*y(0) + \J2 V (- D ^y ° fl )(0) 

= 5 E VIP**®) + 5 E V(-£>_ w y(-e)). 

The last equality follows since for p &1Z we have 

D f> yoR($=yoR(t; + p)- y oR(t;) 

= y(-e - p) - v(-£) 

Employing periodicity of y, we can shift the second summation over 
— JIat back to Qn and upon relabelling, obtain 

**{v) = I E ^(^(0) + \ E n--D-^(0) 
= £ v-(Aiy(0)- 

Finally, point symmetry of TZ and of V are obvious. □ 

Point symmetry of V implies a symmetry of the partial derivatives 
of V in homogeneous states. For p G TZ and g G (IR^) 7 ^, we define 
Vp(gf) := £p-(<7). If V is point symmetric, then we have 

V p (F-n) = -V- p (F-n) for all F G GL(d). (17) 



To see this, one differentiates the identity (15) with respect to g p and 
then evaluates it at g — F ■ TZ. 

We can now prove the that the Cauchy-Born approximation for Bra- 
vais lattices is second order accurate, which was previously observed 
in E, 10, 14 . We nevertheless give a complete proof of the result, since 
it motivates our subsequent analysis for 2-lattices. For Bravais lattices, 
there is no reason to make a distinction between the atomistic and con- 
tinuum variables since we will not consider multiple ways of relating 
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the two. Thus, given a fixed 1-periodic deformation Y : 
let 



-» 



we 



Y N (x) 



NY 



and we interpret Y N as both an atomistic and a continuum deforma- 
tion. 

Proposition 3.2 (Convergence for Bravais lattices.). LetY E C 3 {R d ;W l ) 
be a 1-periodic deformation, and let Y N be defined by (18); then, 

\N~ d g*(Y N ) - £ C (Y)\ < CN' 2 { || V 2 Y ||^ + IIV^IU}, 

where the constant C is a function oflZ and ||V||c 2 - 



Proof. Since the midpoint rule is exact for affine functions, a standard 
quadrature estimate gives 

N~ d S c {Y N ) = N- d f V(V n Y N ) dx 



N- d J2 V{V n Y N (0) 
+ 0(||W(V^)||J 



(19) 



Set V p {g) := £-V{g), and let V p {£) := V p (V n Y N (Z)). By Taylor's 
theorem, 

pen 



+ o(\\v\\ c ,d p y i \o-^ p y n (o n, 



and 



d p y«(z) - v,y A (fl = y»(z + p) - r v (£) - v p f jv (0 

= ^v^ iV (0 + o(||v 3 r Af || ). 

o P v ^' \M Moo/ 



Combining (|20|) and (|21|), we have 

1 
2 



^(A^(0) - v(v*y") = ^W • v ? yiV (0 



pen 



(20) 



(21) 



(22) 



+ o({||v 2 r Af ||L + ||v 3 r iV || 00 }) 



By symmetry of the partial derivatives of V (17), the first term on the 
right hand side of (22) vanishes. We have 



pen pell 

= IJ2{ v p -v 2 p yN -v p -K yN )=u 



pell 
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Figure 1. Point reflection of the honeycomb lattice 
about a lattice site yields a shifted honeycomb lattice. 

where all functions above are evaluated at £ G Qn- Therefore, combin- 

| V 3 F V| 



ing (|19j),(|20j), and (|22j) gives 

N- d \SHY N )-S c (Y N )\ <C{\\V 2 Y N \\ 2 \\V Z Y N \\ ]. 

I V / V / I — U| lloo II llooJ' 

where the constant C depends on 1Z and ||V||c2. Finally, we observe 
that 

|2 



Iv 2 ^!! 2 



iv- 2 ||v 2 r| 



and 



3vJV| 



VY 



nvT 



and the result follows. 



D 



4. Symmetries of the 2-lattice site energy 



We showed in the previous section how the point symmetry (15) 



of Bravais lattices implies second order accuracy of the Cauchy-Born 
energy in the continuum limit. At first glance, 2-lattices do not possess 
a point symmetry: e.g., inversion of the honeycomb lattice about a 
lattice point yields a shifted honeycomb lattice; see Figure [T] Note, 
however, that inverting the lattice about the center of an edge leaves 
it invariant. 

The purpose of this section is to explore how this and other sym- 



metries may be exploited to extend (15). We will identify two such 



non-trivial extensions: for the case of a single species 2-lattice with a 
general site energy, and for the case of a multi-lattices with pair inter- 
action energy (however, we will only formulate the result for 2-lattices). 
In Sections [5] and [6j we will then exploit these observations to derive 
"simple" second order continuum models. 

4.1. Symmetries for single-species 2-lattices. Any single-species 
2-lattice has a reflection symmetry about the "centroid" of a lattice 
site; see Figure § Algebraically, let B e GL(d), peR d \ BZ d , and let 
Jzf := BZ d (J (BZ d + p) be a single-species 2-lattice. Equivalently, we 
may shift the lattice by — p/2 to redefine it as 

JSf := {BZ d -D |J {BZ d + l), 

which immediately reveals the symmetry 

— JSf ' = JSf (23) 
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Figure 2. Point reflection of a 2-lattice about the 
"centroid" x of a lattice site, containing {°,«}. The lat- 
tice remains invariant, however the atom indices are re- 
versed. If the atoms are of the same species, then the 
physical configuration remains invariant. 

Note, however, that if .if has two species, then this operation reverses 
the species and is therefore not a symmetry of a 2-species 2-lattice. 

The symmetry ( [l~5| ) can also be thought of as a permutation of the 
reference list: R(p;a) := (—p;->a) where -iO := 1 and -il := 0, then 
RA = A. We can also define an analogous operation on interactions 
-. : Z d x {0, 1} x {0, 1} ->■ 7L d x {0, 1} x {0, 1}, 

-.fo a, £):=(-# -.a, t9). (24) 

The following Proposition provides a 2-lattice analogy to the observa- 
tion that D p y F = —D_ p y F in Bravais lattices, which was a key step 
to obtain second order accuracy of the Bravais lattice Cauchy-Born 
energy. 

Proposition 4.1. Let F E GL(d) and p = (j>o,Pi) € M 2d , then 

D p y F * = -D^ p y F > p Vp G K. 

Proof. Let y = y F,p and p = (p; a, (3) e TZ, then 

-D^y(0)=-(y^(-p)-y^(0)) 

= -(F(-p)+p^) + (p^ a ) 

= Fp + p^ a - p^ p 

One readily checks by case distinction that p^ a —p^p = pp — p a , which 
concludes the proof. □ 

Definition 4.1 (Point symmetry for a single-species 2-lattice). We 
say that the site potential V is point symmetric (with respect to ->) if 
1Z = ->1Z and 

V({g P } P en) =V(- {g^} p en) V{g p } pen E (R d ) n - (25) 

We show next that, under physically realistic assumptions (invari- 
ance of the energy under permutations and isometries) it is always 
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possible to define a point symmetric site potential for a single-species 
2-lattice. 

Proposition 4.2. Suppose that <§ A given by (pi) is invariant under the 
permutation R of the index list, and under the isometry x h-> —x: 

g*{y) = £*(yoR) = g*(-yoR); (26) 

then, 

^(y)= E v( D *v(Z))> 

where TZ = TZU -^TZ, and V : iW 1 ) 11 -> % 

V({g P } P en) ■= \V{{g P } P en) + \V{{-g^} p ^ n ), 



is point symmetric ( 25 ) 



Proof. Fix y e £2^ and let 7?., V be defined as above. We first show 
that y' := —yoR £ ^y. To that end we must show that the associated 
displacement v! := y' — y B, ° is iV-periodic. To see this, let u := y — y B, °, 
then u is A^-periodic and we have 

<(0 = ~(y o R) a (Z) -B£ = -y^(-0 - B^ 

= -w^(-0 + ££-£?£ = -«^(-0- 

Since U-, a is A^-periodic it follows that u' is A^-periodic. 



The rest of the proof is analogous to the proof of Proposition 3.1 



Employing (26), we have 



S\y) = \S\y) + \S\-y o R) 

= \ E ^(^ (0) + ^ E n-^^(-e))- 

The last equality follows since for p = (p; a,j3)ElZ we have 

D p yoR(t;)=yoR(Z + p,/3)-yoR(Z J a) 
= v(-Z - P, ~v0) - !/(-£, --a) 
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Since y — y B, ° is periodic, we can shift the second summation over 
-Qn back to Qn and, upon relabelling the summation variable, obtain 



^(y) = \ E nm®) + \ E n-^^e - o) 

\ E ViPv®) + \ E V(- D ^y(0) 



2 /_^i v »v»// 2 

= E ^(A/(0), 

where e = (1, . . . , 1) G R d . 

Finally, point symmetry of TZ and of V are obvious. □ 

4.2. Symmetry of pair interaction energies. Suppose that <£' A is 
given by the pair interaction model pi) for a 2-lattice with two different 
species; that is, the site energy V : (M d )^ — > R is defined by 

v(9) = lJ2M\ D P y(0\)- (27) 

(Recall the convention that p = (p,a,/3).) We will now identify a 
symmetry in this model that will play a role analogous to the operation 
-i in the previous section. 

The key observation in this case is that, physically, we should require 
(pa/3 = <f>/3on which motivates the operation ~ : 7L d x {0, 1} x {0, 1} — > 
Z d x {0, 1} x {0,1} by 

~(jr,a,P):=(-fr,P,a). (28) 

We immediately obtain the following result. 

Proposition 4.3 (Symmetry of pair interaction site potentials). Let 
the site energy V be defined by (27) and suppose that TZ = ~1Z and 
(pa/3 = 0/3a f or M a, ft G {0, 1}; then V satisfies the point symmetry 

V{{g P } P en) = V({-<7~pW) V{^} peR G (R d ) n . (29) 

Proof. Under the stated assumptions, 

V{{-9~p}peii) = E^d -9(-rM\) 
pen 

= E ^«(I^(-p;/3,«)D- 
pe^iz 



Relabelling p' = —p, a' = (3,(3' = a we obtain (29 ). □ 

Remark 4.1. 1. Unlike the symmetry -> of a single-species 2-lattice, ~ 
does not arise from any globally defined permutation of the atoms of 
the material. For this reason, one cannot expect that an arbitrary total 
energy will be symmetric under ~, and therefore one cannot hope to 
symmetrize general site potentials. 
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2. The symmetry ~ can be immediately applied to n-lattice pair 
interactions for n > 2. All our subsequent results generalize as well. 



5. Second order accuracy of the Cauchy-Born Energy 

In this section, we give a new rule relating the atomistic and con- 
tinuum variables, and we show that under this rule, the Cauchy-Born 
energy for single-species 2-lattices is second order accurate. Given 1- 
periodic deformation and shift fields Y and P, we first rescale them to 
atomic units 

x 

We then define a corresponding atomistic deformation y N by 



Y N (x) 



NY 



and P 



Ni 



P 



X \ 

n) 



<-N 



Ni 



Z/o(6:=^ V (0-o^ V (0 and 



yf(fl:=y"(£) + -P"(£). (30) 



When (30) is used to connect the atomistic and continuum fields, we 
interpret Y N {^) as the deformation of the centroid of the atoms at site 



£; indeed, the inverse transformation of (30) is 
1 



Y»®:= jW +2/f(0} and P*® := y?(£) - y»{t). 
As in the Bravais lattice case, we need a symmetry of the partial 



derivatives of V (17) to show second order convergence. The following 



lemma establishes the appropriate generalization. 



Lemma 5.1. Suppose that V is point symmetric (25) with respect to 
the permutation operator ->. Let V p be defined by (|8j) ; and let V p (g) : = 
-jP-V(g). Then, for any deformation and shift fields Y and P, we have 



V p {V n (Y,P)(x)) = -V-, p (V*(y,P)(aO). 



(31) 



Proof. Differentiating (25) and (29) with respect to g p yields 

V p {{g p }peTi) = -K p ({-^} p6 r). 
Evaluating at g p = V P (Y, P) and noting that 



-v. P (y, 


P) = -V_ P Y - (-/? - -.a)P 




= V P Y - ((1 - p) - (1 - 




= V P F + {p- a)P 




= v p (y,p), 


we obtain (31). 





a 



))P 



a 



We are now in a position to prove second order accuracy of the 
Cauchy-Born model. 
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Theorem 5.1 (Second order convergence for single species 2-lattices). 
Let Y G C 3 (M. d ; R d ) be a 1-periodic deformation field and P G C 2 (R d ; R d ) 
a 1-periodic shift field. Let y N be the scaled atomistic configuration de- 
fined by (30). Then 

\N- d £ & {y N ) -£ C (Y,P)\ <CN- 2 { \\V 3 Y\\ OQ + || V 2 P||oo 

+ ||V 2 F|| 2 + HVPII 2 j. 

1 II v lloo ' II v lloo J 



where the constant C is a function of 1Z and ||V| 



c 2 - 



Proof. Since the midpoint rule is exact for affine functions, a standard 
quadrature estimate gives 

N -dgc ( y tf pTV) =N -d f y (y^ (y N , P N ) ) 

= N- d J2 V (Vn(Y N ,P N )) 

+ 0(\\V 2 V(Vn(Y N ,P N ))\\J. (32) 

It is easy to see that 

||vV(v*(y*,p»))L 

<C(||V 3 F 7V || + \\V 2 P N \\ +\\V 2 Y N \\ 2 +\\VP N \\ 2 ). 

— VII lloo II lloo II lloo II lloo/ 

We now estimate the error in V(Vti(Y n , P n ) (£)) . By Taylor's the- 



orem, 



v(D ny N ) - v(v n (Y N , />*)) = £ v i D P y - v »y} 

pen 

+ 0(\\V\\ c i\D p y N -V p (Y N ,P N )\ 2 ), 

(33) 

where all functions above are evaluated £ G VLn, and V p := V p ( V^ (Y N , P N ) (£)) . 
By a second application of Taylor's theorem, 

D P y N (0 - V p (y*, P N ) (0 = tf{t + P ) - y»(£) 

-V p Y N (Z)-((3-a)P N (0 
= Y N (Z + p) + {(3-l)P N (Z + p) 

-y n (0-(*-1)p n (0 

-V p Y N (0-(P-a)P N (0 



p 

2 



+ 0(||V 3 F 7V || +\\V 2 P N \\ ). (34) 

Ml lloo ! lloo/ "> ' 
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Substituting (34) into (33), and assuming again that all functions are 
evaluated at £ G Qn, we obtain 

V(D n y N )-V(V n (Y N ,P N ))=J2v P -<P) 

pen 

+ 0(||V 3 F Ar || + \\V 2 P N \\ ), (35) 

VII lloo II II OO/ ' v ' 

where 

e(p):= l 2 VlY N (Z) + ({3-l)V p P N (0. 
We now observe that 

ehp) = §V i,Y*(0 + (^ - l)V. p P N (0 

= ^ 2 - P Y N (o + (a-^-^- P p N (o 

= e(p), (36) 

from which we deduce that the the first term on the right-hand side 
of (35) vanishes: using (37) and point symmetry of the derivatives (31 ), 
we obtain 

E v p ■ e ^) = \ E v p ■ e ^) + v *p ■ <*p^ 



pell 



pell 



\Y, V p-{<p)-<*p)) 



pell 



(37) 



Finally, combining (|32j), (|33j), (J34J), (|35j), and (|37j) gives 
\N- d £ a (y N ) -£ C (Y,P)\ 

<c(\\v 3 y n \\ +||v 2 p Ar || +||v 2 r JV || 2 4 

— Ml lloo II lloo II lloo 

The result follows by rescaling 

Hv 3 ^!! =iv- 2 |lv 3 >1l , 



VP 



N\\2 ^ 



l V 2yiV||2 



\V 2 P N \ 



N~ 2 \\V 2 P\ 



and 



VP 



oo 

AMI 2 
oo 



N- 2 \\V 2 Y\Y 

II lie 

iv- 2 llvp|| 2 



D 



Remark 5.1. Suppose that we are in the multi-species 2-lattice case, 
and must use the symmetry ~ instead of ->. In this case, we have 

e(~p) = ^- P Y N (0 + (a - l)V. p P N (0 

= \V 2 p Y N {i) + {\-a)V p P N {Z)^e{p), 

whenever a = 0. Thus, in our above proof the first order terms do not 
cancel and the Cauchy-Born model is only first order accurate in this 
case. We will show in the next section how a new strain energy density 
can be constructed that is second order accurate in this case. 
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6. A SECOND ORDER ACCURATE STORED ENERGY DENSITY FOR 

PAIR INTERACTIONS 

In this section, we show how the symmetry ~ for multi-lattice pair 
interactions can be exploited to derive a "simple" second order accurate 
stored energy density, which does not coincide with the Cauchy-Born 
strain energy density To maintain a notation consistent with the rest 
of the paper, we state all results only for 2-lattices, but stress that 
they can be immediately extended to general multi-lattices composed 
of more than two Bravais lattices. 

Suppose we are again given 1-periodic deformation and shift fields 
(Y, P). Let (Y N , P n ) and y N = (y^ , y±) be defined as in the classical 



Cauchy-Born setting in (11) 



To exploit the symmetry ~ defined in < \28\j , we define a different 
multi-lattice directional derivative. If, for the moment, we interpret 
y N as a continuum field, then we may define 

Vpi^^lVpi^ + IVp^ + d^-rf). 

We observe that V^ p = — V p for this definition. Written in terms of 
(Y N ,P N ), V p becomes 

V P {Y N , P N ) := V p Y N + {p - a)P N + ^fV p P N . (38) 

Note that, in contrast to the previous multi-lattice directional deriv- 
ative (J8J), this is not a scale-invariant field. Indeed, defining, e := l/N 
and rewriting the directional derivative in terms of (Y, P), we arrive at 

V p (y, P) := V P (Y N , P N ) = V P Y + (/3 - a)P + ^eVpP, 

which gives rise to the new stored energy density 

W e (VY,P):=V(V £ n (Y,P))- 

Remark 6.1. The fact that W £ employs shift gradients but not strain 
gradients makes the resulting model straightforward to implement us- 
ing C°-conforming discretizations. By contrast, the higher-order mod- 
els derived, e.g., in fl] require also higher-order conforming discretiza- 
tions. 

The symmetry V^ p = — Vp was the first key ingredient in the proof 
of of second order consistency of the single-species 2-lattice strain en- 
ergy density. We now generalize the second key ingredient, the sym- 
metry of partial derivatives of V . 



Lemma 6.1. Suppose that V is point symmetric (29) with respect to 
the permutation operator ~. Let V p be defined by (|8j), and let V p (g) : = 
i^-V(g). Then, for any deformation and shift fields Y and P , we have 

V p (V n (Y N ,P N )(x)) = -V^(V n (Y N ,P N )(x)). (39) 



Proof. The proof closely follows the proof of Lemma 5.1 □ 
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Theorem 6.1 (Second order convergence for 2- lattices). Let Y G 
C 3 (M d ;]R d ) be a 1-periodic deformation field and P G C 2 (M d ;IR d ) a 1- 
periodic shift field. Let y be the scaled atomistic configuration defined 



oo 



by (11); then 

N- d \^{y N ) -£ C (Y N ,P N )\ < CN~ 2 { llV^Hoo + ||V 3 P| 

+ iiv 2 pii ■+- nv 2 yn 2 + iivpii 2 \ 

T || V -T ||oo "I" || V I || 00 T^||V-T|| 00 J'. 

where the constant C is a function oflZ and ||V||c 2 - 



Proof. The proof of the result is the same as the proof of Theorem 5.1 



up to (33), except that we need P G C 3 for the quadrature estimate, 



since VP now enters the definition of the multi-lattice directional de- 



rivative. A computation analogous to (34) yields 

D p y N - V P (Y N , P N ) = ^V P P N + \V 2 p Y N 

+ o{\\v 3 y n \\ 00 + \\v 2 p n \\ 00 ). 



Continuing as in the proof of Theorem |5.1[ we obtain 

V(D n y N )-V(Vn(Y N ,P N ))=J2 V f<P) 

pen 

+ 0(||V 3 Y wv || +\\V 2 P N \\ ), (40) 

V II lloo II II oo/ ' v / 

where 

e(p) := ^V P P N + \V 2 p Y N . 

It is straightforward to see that e(~p) = e(p), and hence the rest of 



the proof follows Theorem 5.1 □ 



7. Conclusion 

We have identified two new second order continuum models for multi- 
lattices. Our approach is based on the identification of symmetries 
similar to the point symmetry of Bravais lattices. We then extend the 
standard proof of second order accuracy of the Cauchy-Born rule for 
Bravais lattices to derive second order models for multi-lattices. 

For single-species 2-lattices, we show that the classical Cauchy-Born 
model is of second order in the continuum limit, provided that the 
atomistic and continuum kinematic variables are related in a new way. 
If the classical relationship between the variables is adopted, then the 
Cauchy-Born model is only first order. We also give a new stored en- 
ergy density for a general multi-lattice modeled using pair interactions, 
and we show that this energy is second order accurate, as well. 

These results are being used to develop accurate atomistic/continuum 
couplings. We show in [22] that the interfacial error of a coupling can 
be dramatically reduced when a point symmetric site energy is used, 



and we make similar applications of the results of this paper in 17 . We 



also remark that our methods achieve second order accuracy without 
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requiring the use of ^-conforming elements, and so our models are 
easier to implement than the second order models of fl]. 
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